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Abstract 

A  simple  two-level  radiation-scaling  model  has  previ¬ 
ously  been  shown  to  agree  with  0.7-MA  neon  gas-puff 
and  7-MA  aluminum  wire-array  experiments.  Here,  the 
scalings  of  optimum  load  mass  and  peak  K-shell  yield 
with  current  for  argon  gas  puffs  are  compared  with  data 
over  an  order-of-magnitude  current  variation.  Optimum 
load  conditions  are  determined  by  maximizing  K-shell 
yield  subject  to  fixed  implosion  energy,  resulting  in  peak 
yield  and  optimum  mass  as  functions  of  implosion  energy 
and  stagnation  radius.  The  scalings  with  current  result 
from  implosion-energy  and  implosion-time  relations  that 
fit  Double  Eagle  nested-shell  argon  gas-puff  results.  Both 
the  maximum-yield  and  optimum-mass  scalings  provide 
good  absolute-value  fits  to  the  data,  including  the 
variation  with  implosion  time. 

I.  BACKGROUND  AND  MODEL  REVIEW 


tion  radius.  The  current  scalings  result  from  the  im¬ 
plosion  energy  vs  current,  implosion-time  scaling,  and 
compression  ratio  that  fit  Double  Eagle  (DE)  nested-shell 
argon  results.  [5] 

The  two-level  model  is  a  simple,  first-principle  treat¬ 
ment  useful  for  estimating  how  K-shell  x-radiation  from 
the  PRS  scales  with  plasma,  pulsed-power,  and  imploding- 
load  parameters.  [1]  Input  parameters  for  the  model  are 
the  implosion  energy  Ew(J/cm),  load  mass  m {g/cm),  stag¬ 
nation  radius  Rf(cm)  of  the  radiating  plasma,  and  atomic 
number  ZA.  The  only  radiation  considered  is  the  line  asso¬ 
ciated  with  the  principal  2-to-l  transition  of  a  single,  K- 
shell  ionization  state,  calculated  for  uniform  density  and 
stagnated-plasma  temperature  T(eV)  within  Rf.  The  tem¬ 
perature  derives  from  a  simple  energy-balance  equation. 
For  gas-puff  z-pinches,  the  implosion  energy  is  calculated 
from  local  snow-plow  implosions  in  measured  gas-density 
distributions  and  contains  both  radial  kinetic  energy  and 
internal  energy  due  to  shock  heating.  [6] 

Energy  balance  takes  the  form  [  1  ] 


A  simple  radiation- scaling  model  based  on  two-level 
atomic  physics  and  plasma  energy  balance  [1]  was  de¬ 
veloped  to  guide  selection  of  plasma-radiation-source 
(PRS)  load  parameters  and  estimate  K-shell  x-ray  yields 
for  z-pinch  drivers.  The  model  has  been  shown  to  agree 
in  detail  with  neon  gas-puff  experiments  [2]  at  0.7  MA 
and  aluminum  wire-array  experiments  [3]  at  7  MA.  Here, 
it  is  used  to  predict  the  current  scaling  of  peak  K-shell 
yield  from  optimum  load  parameters  for  argon  gas-puffs. 
Results  are  compared  with  data  from  a  variety  of 
generators  and  load  conditions  spanning  an  order-of- 
magnitude  current  variation,  including  recent  results  on 
the  Sandia  National  Laboratories  Z  generator  at  15  MA 
[4].  In  light  of  the  variety  of  driver  and  load  configura¬ 
tions  represented  by  the  data,  it  is  noteworthy  that  both 
the  maximum-yield  and  optimum-mass  scalings  provide 
good  absolute-value  fits  to  the  data,  including  the  implo¬ 
sion-time  variation. 

Optimum  temperature  and  mass  are  determined  by 
maximizing  K-shell  yield  subject  to  energy  balance  with 
fixed  implosion  energy,  resulting  in  peak  yield  and  opti¬ 
mum  mass  as  functions  of  implosion  energy  and  stagna¬ 


EW=EI+Y  +  Ep  (1) 

with  all  quantities  in  J/crn  of  load  length.  In  Eq.  (1),  Ew  is 
the  implosion  energy,  E:  is  the  stagnated-plasma  internal 
energy,  Y  is  the  K-shell  radiated  energy,  and  Ep  contains 
energy  losses  in  neglected  processes  such  as  radiation  from 
other  excitations  and  residual  energy  in  the  plasma 
rebound  following  peak  compression.  The  internal  energy 
in  the  K-shell  regime  can  be  written 

E,  =1.8xl05mT/Z“s  (2) 


The  radiation  term  in  Eq.  ( 1 )  is  given  by  Y  =  S Yt ,  where 
the  optically-thin  radiation  rate  is  given  by 


Yt  = 


3.9xl016m2 

RfTZA2 


-exp 


(  10.2Z2  A 


(3) 


and  S  is  the  radiation  escape  fraction.  For  the  conditions 
of  this  study,  S  =  1  for  plasma  parameters  that  maximize 
the  yield.  [1]  The  model  is  in  good  agreement  with  neon 
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and  aluminum  experimental  yields  when  Ep  is  taken  to  be 
Ew/2,  so  that  energy  balance  becomes 

E\v  /2  =  E,  +  Yt  (4) 

It  is  convenient  to  normalize  the  load  mass  and  stagna¬ 
tion  temperature  according  to 

y  =  T/T*  ;  n  =  m/m*  (5) 

where 

T*  =  3.4Z^  (6) 

m*  =  8.3xlO_7Ew  /Z^72  (7) 

so  that  Eqs  (2)  and  (3)  take  the  form 

l  =  ny  +  2Yt/Ew  (8) 


The  variations  of  nop  and  Ymax/Ew  with  Ew/Kbi  are  shown 
in  Fig.  1.  The  asymptotic  behaviors  in  the  weak-emission 
(I4)  and  strong-emission  (I2)  regimes  are  indicated. 
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(9) 


where 

Kb  (J/cm)  =  Rf Z^64  /  390  (10) 

characterizes  the  transition  energy  centered  between  the 
asymptotic  I4  and  I2  regimes.  [1] 


Figure  1.  Variations  of  normalized  load  mass  and  K-shell 
efficiency  with  implosion  energy. 

The  model  variables  are  next  related  to  implosion  time 
tjmp  and  Ew  in  order  to  recover  the  scalings  with  current 
from  the  general  results  of  Eqs.  (11)  -  (13).  Electrical- 
energy  balance  and  the  equation  of  motion  provides  the 
following.  [1] 


II.  YIELD  OPTIMIZATION 


The  objective  is  to  maximize  Y,  constrained  by  energy 
balance  for  fixed  implosion  energy  Ew  .  The  implosion 
energy  depends  on  the  load-current  history  (and  therefore 
the  implosion  time)  and  the  ratio  R0/Rf  ,  where  R0  is  a 
measure  of  the  outer  radius  of  the  gas  distribution.  The 
implosion  time  is  proportional  to  mR0  ,  so  that  constant 
Ew  requires  R0  ~  m'1/2,  and  fixed  compression  ratio  re¬ 
quires  Rf  ~  R0 .  For  fixed  Ew ,  the  radiation  term  takes  the 
form 


E  n5/2 

n 


Kk 


-exp 


3  2 

y. 


(11) 


where  Kbl  is  for  Rt  =  Rn  corresponding  to  n  =  1. 

Maximizing  Eq.  (11)  subject  to  Eq.  (8)  results  in  the 
optimal  stagnation  temperature  given  by  yop  =  6/7  and 
load  mass  nop  given  by 


0-7l(l-  6n  !l) 


512 
1  op 


(12) 


for  which  the  maximum  K-shell  yield  is 


/  D  \ 

Ew(7/c/w)  =  103yIo(M4)  •  In  — ^ 

lR‘  j 

-6 mU2(g/cm)R0(cm) 


timpC-5)  =  10  Ct 


I  o(MA) 


(14) 

(15) 


The  three  constants  y=0J,  Ro/Rf  =  14,  and  Ct  =  17  are 
chosen  to  fit  DE  argon  nested-shell  experiments,  where  K- 
shell  yield  vs  implosion-time  was  obtained  by  varying  the 
nozzle  pressure  (load  mass)  [5].  The  data  are  seen  in  Fig. 
2  to  compare  well  with  the  model  when  the  mass  and  peak 
current  of  each  shot  were  used  to  compute  individual  mo¬ 
del  yields.  Two  radiation  calculations  are  performed:  one 
using  the  full  implosion  energy  (kinetic  +  shock  heating) 
that  agrees  with  the  data,  and  a  second  using  only  the  ki¬ 
netic  portion  that  gives  about  half  the  measured  yields. 
The  results  demonstrate  that  shock  heating  should  be  in¬ 
cluded  in  radiation  modeling  to  predict  the  performance  of 
a  gas-puff  or  nested- wire  PRS. 

Substituting  the  three  constants  into  Eqs.  (14)  and  (15) 
leads  to 

Ew  (J/cm)  =  1 850  •  Iq  ( MA)  (16) 


R0op(cm)  _  1.8X107 

W*)  “  <2 


(17) 


|(l-6nop/7) 


An  outer  radius,  determined  from  Eq.  (17)  with  n  =  1,  is 
(13)  used  to  determine  Rfl  from  the  compression  ratio,  which  is 
then  used  to  calculate  Kbl  from  Eq.  (10).  The  variation  of 


peak  argon  K-shell  yield  with  current  and  implosion  time, 
is  then  determined  from  Eqs.  (12),  (13)  and  (16). 
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Figure  2.  K-shell  yield  vs  implosion  time  for  DE  argon 
nested-shell  loads  compared  with  the  model. 
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III.  RESULTS 

The  variation  of  predicted  maximum  argon  K-shell 
yield  with  current  is  shown  in  Fig.  3  for  100-ns  and  300- 
ns  implosion  times.  The  yield  smoothly  varies  from  the 
expected  I4  to  l2  dependence  as  the  current  increases,  with 
the  center  of  the  transition  at  about  8  MA  for  argon.  For 
other  atomic  numbers,  the  definition  of  Kb  provides  a 
transition  current  that  varies  like  ZA3'3.  In  the  I4  regime. 


Figure  3.  Maximum  argon  K-shell  yield  vs  peak  current 
for  various  generators  and  loads  compared  to 
the  model  for  timp  =  100  and  300  ns. 


the  model  yield  varies  inversely  with  implosion  time.  At 
higher  currents,  this  long-implosion-time  yield  penalty  is 
reduced,  suggesting  that  high-current,  long-implosion-time 
drivers  may  achieve  yields  approaching  those  of  100-ns 
drivers. 

Also  plotted  in  Fig.  3  are  the  maximum  argon  K-shell 
yields  achieved  with  various  generators  and  loads.  The 
labels  “short”  and  “long”  correspond  respectively  to  im¬ 
plosion  times  in  the  70-  to  110-ns  and  190-  to  300-ns  re¬ 
gimes.  Circles  and  triangles  correspond  respectively  to 
uniform-fill  and  nested-shell  nozzles. 

The  variation  of  predicted  optimum  argon  mass  with 
current  is  compared  to  various  measurements  in  Fig.  4. 


Figure  4.  Optimum  argon  load  mass  vs  peak  current  for 
various  loads  compared  to  the  model. 

The  circles  are  kinetic  masses  based  on  matching  snow¬ 
plow  implosions  to  observed  implosion  times.  The  error 
bar,  masses  determined  by  interferometry  for  the  DE 
nested-shell  experiments  [6],  agree  with  the  kinetic  values. 
The  square  is  an  estimate  of  mass  for  the  highest  K-shell 
yield  Z  shot.  [4]  The  data  show  that  the  optimum  mass 
does  not  depend  on  implosion  time,  in  agreement  with 
modeling.  The  dashed  line  extends  the  weak-emission 
scaling  to  higher  current.  The  optimum  mass  is  propor¬ 
tional  to  I2  for  low  currents  and  transitions  to  I12,  in 
agreement  with  the  higher-current  data. 

The  weaker  current  scaling  of  mass  for  strong  radiators 
should  be  considered  in  the  design  of  high-current  loads, 
though  the  predicted  yield  reduction  for  non-optimal  load 
masses  is  weaker  than  for  low  currents.  If  so,  the  mass  at 
high  current  may  be  increased  (initial  radius  decreased)  for 
mediation  of  the  Raleigh-Taylor  (RT)  instability  with  a 
minor  K-shell  yield  penalty.  This  possibility  is  illustrated 
in  Fig.  5,  where  the  relative  yield  vs.  load  mass  is  plotted 
for  argon  at  low  and  high  current,  corresponding  to  Ew/Kbi 
=  0.1  and  3.6.  At  high  current,  load  masses  a  factor-of-two 
greater  than  optimum  are  predicted  to  reduce  K-shell  yield 


by  about  10%.  For  aluminum  loads,  the  corres-ponding 
low  and  high  currents  in  Fig.  5  are  reduced  to  0.85  and 
5.1  MA,  and  the  weak  high-current  dependence  of  yield 
on  masses  above  optimum  has  been  demonstrated  [1,  3]. 
Based  on  a  single  non-optimum  argon  shot.  Fig.  5  shows 
that  this  trend  has  not  yet  been  observed  on  Z.  [4] 

The  possibility  of  using  masses  more  massive  than 
optimum  with  minor  yield  penalty  may  be  most  useful  for 
high-current,  long-implosion-time  drivers  such  as  the 
Decade  Quad  and  Dual  Quad,  where  timp  =  300  ns  for 
argon  leads  to  an  optimum  load-radius  of  about  6  cm. 
Load-performance  risk  reduction  may  be  achieved  by 
decreasing  R0opt  to  about  4  cm  by  a  factor-of-two  increase 
in  load  mass  above  optimum. 


m/niopt 

Figure  5.  Variation  of  argon  K-shell  yield  with  load 
mass  for  currents  in  the  weak-  and  strong- 
emission  regimes. 


IV.  CONCLUSIONS 

A  simple  radiation-scaling  model  based  on  two-level 
atomic  physics  and  plasma  energy  balance  has  been  used 
to  determine  the  current  scalings  of  optimum  load  mass 
and  peak  K-shell  yield  for  Ar  gas  puffs.  Results  of  the 
model  were  compared  with  data  over  an  order-of-magni- 
tude  current  variation.  In  light  of  the  wide  variety  of 
driver  and  load  configurations  sampled,  it  is  noteworthy 
that  both  the  maximum-yield  and  optimum-mass  scalings 
provide  good  absolute-value  fits  to  the  data,  including  the 
implosion-time  dependence.  Agreement  with  the  data  de¬ 
pends  on  the  use  of  the  full  implosion  energy  (kinetic  + 
shock-heated  internal)  of  distributed-mass  loads. 

In  addition  to  extending  validity  to  atomic  number  18 
from  previous  neon  and  aluminum  benchmarks,  agree¬ 
ment  of  argon  data  with  the  zero-dimensional  radiation¬ 


scaling  model  also  suggests  that  gas  puff  performance  is 
not  greatly  degraded  by  instability  effects.  This  obser¬ 
vation  supports  the  snowplow  mechanism  proposed  for 
suppression  of  RT  modes  in  distributed-mass  loads.  [7] 
Work  supported  by  DTRA. 
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